Two-scale localization in disordered wires in a magnetic field. 
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Calculating the density-density correlation function for disordered wires, we study localization 
properties of wave functions in a magnetic field. The supersymmetry technique combined with the 
transfer matrix method is used. It is demonstrated that at arbitrarily weak magnetic field the far 
tail of the wave functions decays with the length Leu = 2Lco, where Leo and Leu are the localization 
lengths in the absence of a magnetic field and in a strong magnetic field, respectively. At shorter 
distances, the decay of the wave functions is characterized by the length L^o- Increasing the magnetic 
field broadens the region of the decay with the length Leu, leading finally to the decay with Leu 
at all distances. In other words, the crossover between the orthogonal and unitary ensembles in 
disordered wires is characterized by two localization lengths. This peculiar behavior must result 
in two different temperature regimes in the hopping conductivity with the boundary between them 
depending on the magnetic field. 
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Disordered systems have been under intensive study for 
several decades but only recently strong localization at 
moderate disorder in ID wires has become the subject of 
a systematic experimental investigation 0| . The authors 
of Ref. studied electron transport in submicrometer- 
wide wires fabricated from Si (5-doped GaAs. A large 
number of the wires were connected in parallel and the 
conductivity of the entire system was measured as a func- 
tion of temperature T and magnetic field H applied per- 
pendicular to the wires. The activation temperature de- 
pendence observed allowed to demonstrate the exponen- 
tial localization and extract the dependence of the local- 
ization length Lc on the magnetic field. It was found that 
in a strong magnetic field the localization length is twice 
as large as in the absence of the magnetic field. 

Localization due to quantum interference is a funda- 
mental property of disordered one-dimensional (ID) sys- 
tems. Mott and Twose predicted that at arbitrary 
weak disorder all states of a disordered chain had to be- 
come localized, which has been proven later |^ . Thouless 
Q realized that all wave functions had to be localized also 
in thick wires. Such wires are more interesting because 
electrons can diffuse at distances exceeding the mean 
free path I but get localized on the localization length 
Lc = I {pqS) , where po is the Fermi-momentum and S is 
the cross-section of the wire. Explicit calculations con- 
firmed this picture for thick wires ||^,^ as well as for sys- 
tems of coupled chains 0. In contrast to the chains, the 
electron motion in thick wires is sensitive to an external 
magnetic field. A remarkable effect, namely, doubling of 
the localization length when applying a strong magnetic 
field H was predicted for such systems 7|. It is this 
doubling that has been observed in Ref. [1|. Different 
aspects of the localization in wires as well as connections 
to random matrix theories were discussed later [p|-p^ . 



Surprisingly, little attention has been paid to the 
crossover between the limits of zero and strong magnetic 
fields H (between the orthogonal and unitary ensembles). 
Study of this crossover is important, however, because 
the localization length can be measured experimentally. 
Except for an interpolation formula for the localization 
length suggested in Ref. ||l^ and a numerical study of 
Ref. practically no attempt has been made to de- 
scribe the crossover. Apparently, the absence of atten- 
tion has been due to a common belief that the crossover is 
simple and the only thing that can happen is that the lo- 
calization length changes smoothly between its value Leo 
at zero field and Leu = 2Lco at a strong magnetic field. 
Such a scenario was used, for example, by the authors of 
Ref. ^ in their attempt to fit the data on the interpo- 
lation curve of Ref. |13|, suggested basically to describe 
the smooth change of the localization length. 

In this Letter we present results of an analytical study 
of a correlation function describing the decay of wave 
functions in disordered quantum wires. In contrast to the 
one-scale picture, the behavior of the wave functions at 
finite magnetic fields turns out to be more complicated. 
The most striking picture of this decay is that, even at 
very weak magnetic fields, the far tail of the wave func- 
tions falls off with the length Leu, whereas at smaller 
distances the decay is governed by the length Leo. The 
larger length Leu does not depend on the strength of the 
magnetic field. Increasing the magnetic field results in a 
broadening of the asymptotic region until the entire wave 
function starts decaying with the length Leu- 

The behavior of localized wave functions can be well 
described by the correlation function poo (f) 
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where ipa {t) and Eq are the eigenfunction and eigenen- 
ergy of a state a, r > is the coordinate along the wire, 
and the angle brackets stay for averaging over disorder. 

The function poo {f) is important not only for a the- 
oretical description of localized wave functions but it 
also determines directly the hopping conductivity at low 
temperature. According to results obtained quite long 
ago , the hopping conductivity of one-dimensional 

chains and wires is described by a simple formula 



vi^vS, (2) 



where v is the density of states, vi is the one-dimensional 
density of states, S is the cross-section of the wire and 
Lc is a localization length that can be extracted from the 
function (r) , Eq. ([I]) . The pre- factor ctq depends on 
parameters of the model. Remarkably, even the Coulomb 
interaction does not change the temperature dependence, 
Eq. (P), entering ctq only |15|. Equation (0) was used in 



Ref. |1| for extracting the localization length from the 
temperature dependence of the conductivity. 

The correlation function p^o {t) can be found using the 
supersymmetry technique P,pT|. Let us first discuss the 
final result in a weak field which reads 
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where Leo — t^i^iDq is the localization length for the or- 
thogonal ensemble, icu = 2ico is the localization length 
for the unitary ensemble, Dq is the classical diffusion co- 
efficient, and a is a positive constant of order unity. The 
parameter X = 2'K(j)/(j)Q describes the crossover between 
the orthogonal and unitary ensemble, 0o = hc/e is the 

1/2 

flux quantum, and cf) — 2HLco (y^) is a characteristic 
magnetic flux through an area limited by the localization 
length. The coordinate y is perpendicular to the direction 
of the wire and the magnetic field. The symbol 
stands for the averaging across the wire. For the "fiat" 
wires made on the basis of a 2D gas as in Ref. |l| , one has 
(y^) = d/\/T2, where d is the width of the wire. For 

1 /2 

wires with a circular cross-section {y'^') — d/A, where 
d is the diameter. Equation (|^) is written in the limit 
X < 1 and r > Leo. At X = 0, Eq. (|) reduces to 
the well known result for the wire without a magnetic 
field PJO]. At small but finite X <C 1, the second term 
in Eq. (0) is small for not very large r but in the limit 
r — > oo it is always larger than the first one. 

Equation (|^) shows that a weak magnetic field changes 
drastically the tail of wave functions leaving their main 
body almost unchanged. Even if there is a correction to 
the localization length Leo, the second term in Eq. (|^) is 
more important in the limit of large r. In other words, 
the characteristic behavior of wave functions at distances 



r <C rx is described by formulae of the orthogonal en- 
semble, while the tail r ^ rx corresponds to the unitary 
one. Comparing the first and the second term in Eq. ^ 
one estimates the characteristic distance rx 



rx 
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The fact that any weak magnetic field changes the tail 
of the wave function looks quite natural because regions 
where the amplitude of the wave function is very small 
must be very sensitive to external perturbations. The 
localization of the wave functions is a result of multi- 
ple interference of waves scattered by impurities. The 
tails are formed by a coherent propagation of waves along 
very long paths, which makes the influence of any weak 
magnetic field on the interference extremely important. 
The separation into the "orthogonal" and "unitary" parts 
loses its sense at X ~ 1 corresponding to characteristic 
magnetic fields determining the crossover between the 
orthogonal and unitary ensembles. 

Analogous changing of the asymptotic behavior in a 
magnetic field occurs in OD (see, e.g. The level- 

level correlation function for the orthogonal ensemble is 
proportional to the energy u in the limit w — s- but 
any weak magnetic field changes this behavior to an 
-dependence, which is characteristic for the unitary en- 
semble. Of course, the region of the lj^- dependence is 
very narrow at weak magnetic fields. 

The factor X^ In^ X in Eq. (j^) is proportional to the 
square of the flux through a segment of the wire limited 
by rx, Eq. (^) and resembles the Mott law for 
the ac conductivity, (t{u)) ^ uP' \t? lo. Using the standard 
perturbation theory one can understand that rx is the 
length at which the magnetic field considerably changes 
correlations and so, its presence in Eq. is natural. 

The final result, Eq. (^, is essentially non-perturbative 
and its derivation is not simple. Let us present now the 
scheme of the computation. 

The correlation function poo {t) is extracted from the 
low-frequency limit of the density-density correlation 
function which, in its turn, is presented as a functional 
integral over 8x8 supermatrices Q [O 
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(r) = -27r2j.2 J (Q) Qj] (r) exp {-F [Q]) DQ 

where Qii ^^'^ Qil ^-re the matrix elements of superma- 
trices Q. The free energy functional F [Q] is written in 
ID in the presence of a magnetic field as 



F[Q] = ^Str 



DoidQ{x)f + 2iujAQ{x) dx . (6) 



In equation (j6|), we introduced the operator dQ = VxQ — 
{ie/hc)A[Q, Ta], the coordinate x is chosen along the wire, 
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and the standard notations for the supertrace Str and 
matrices T3 and A are used |11|. Choosing the gauge 
A = {Hy, 0, 0) and using the notations introduced below 
Eq. (^) one rewrites the functional F as 

F[Q]^F [Q]^^o + (32Lco)"' Str J[Q, t^]^ dx . (7) 

Calculation of the one-dimensional functional integral 
in Eq. (|^) is performed using the transfer matrix tech- 
nique. The Fourier transform of the correlator is writ- 
ten as the integral 0| 



(k) = 



(8) 



of the product of the scalar and the matrix functions 
'^{Q) and P{Q), representing the partition functions of 
the segments a; < and x > 0, respectively. Due to the 
supersymmetry one has / ^(Q)dQ = 1. 

Discretizing the wire and using Eqs. ^ one de- 
rives recurrence relations for the functions 5" and P on 
neighboring sites. In the continuum limit, these relations 
yield effective "Schrodinger equations." The correspond- 
ing Hamiltonian contains the Laplacian Ag acting in the 
space of matrix elements The functions ^' and 

P have the symmetry of the states with the zero and the 
first "angular momentum," respectively. In a schematic 
form, one has 



- c;;Str(Ag)] * = , (9) 
wStr(AQ)] Pfe - Q * , 



where Aqq and Aig are the projections of the Lapla- 
cian Aq on the states with the zero and first angular 
momentum, and uj = iix)L^^/2DQ. 

Equations (|[ ^ solving the problem are rather gen- 
eral. A certain parametrization of Q-matrices should 
be chosen to perform explicit calculations. As soon 
as H 7^ 0, the standard parametrization used in Refs. 
|p|,|6Ulll] is not convenient, so that we use the "magnetic 
parametrization" of Refs. ||l^,|ll|] constructed for study- 
ing the crossover between the orthogonal and unitary en- 
sembles in OD. However, even this parametrization, when 
applied to Eqs. (||) leads to extremely complicated par- 
tial differential equations. 

Instead of trying to solve these equations exactly, we 
concentrate on a study of the limit of weak magnetic 
fields. The localization length is determined by the poles 
of the function P in the fc-plane ||ll[] in the region of 
"free motion" , faraway from the "barrier" given by the 
last term in Eq. (j^). The position of the poles fc„ are 
determined by the equation £n(fc„) — 0, where £'n(fc) are 
eigenvalues of the operator entering the LHS of Eq. (|^) 
corresponding to the eigenfunctions (pn{Q)- According 
to the general procedure developed in Ref. one has 
no need to consider solutions at arbitrary uj. Although 



the dependence on lo is necessary for computing different 
matrix elements, the poles can be found just by putting 
Lu = 0. This is natural because the localization prop- 
erties are determined for finite samples at w = [ p^ . 
Even then, we cannot find all solutions Pk (Q) and cor- 
responding poles k. However, in order to determine the 
function 1^ (r) at large distances, we need only one state 
with the smallest non-zero kg (X). At zero magnetic field, 
such a state ipo (Q) is well known [^,0 to have a pole 
at A;o (0) = ±z (4Lco) This corresponds to the first 
term in Eq. . As soon as one applies a weak magnetic 
field, the state is distorted. This effect is not crucial: 
Wc find no corrections in the lowest order either to the 
localization length or to the pre-exponential. 

The main result of the present work, namely, the sec- 
ond term in Eq. (^, can be obtained without complicated 
calculations. It turns out that, at arbitrarily weak mag- 
netic field, an additional state (pi {Q) with a smaller value 
ki = zti (8Lco)~ appears. It is this state that determines 
the behavior of Poo (f) in the limit r —^ 00. Remarkably, 
the value fci does not depend on the magnetic field even 
in the limit X 00, when the state ipi gives the main 
contribution at all distances. 

The origin of the states, ipo and tpi, is apparent. An es- 
sential feature of the magnetic parametrization is a finite 
contribution arising due to the singularity of the Jaco- 
bian at Aic, Ac ^ 1, where Aic and Ac are "eigenvalues" 
corresponding to the cooperon degrees of freedom. This 
singularity is of the type (Aic — Ac) ^ and is usual in the 
supersymmetry technique. The procedure of regulariza- 
tion of the singularity has been developed for the OD cases 
in Refs. p6 11 1. In the ID case, the contribution from the 



singular term gives the correlator for the unitary ensem- 
ble (in the limit Aic, Ac ^ 1 the cooperon degrees of 
freedom are frozen). In order to compensate this part at 
moderate fields, the regular contribution should contain 
both the solutions (po and (pi. At X ^ 0, the solution 
ifi exactly compensates the singular term, so that the 
orthogonal part is given by ipo only. 

Having disregarded the term with the external fre- 
quency UJ, we arrive at a Hamiltonian with separated 
cooperon and diffuson variables. However, it is still too 
complicated since essential values of Ac and Aic are such 
that Ale, Ac — 1. This difficulty is avoided by considering 
the derivative over the magnetic field, Y^. For the quan- 
tity Y^, in contrast, large Aic and Aid become important, 
which substantially simplifies the "Hamiltonian" 

w = -^IddL - ^Irdl^ + 2AicaA,. + x'xl . 

The solution decaying at Aic — * 00 equals 

4- = (1 XAic) exp(-XAic) , (10) 

and satisfies the boundary condition ^'(0) = 1. The main 
contribution to the correlator Y^ comes from the region 
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1 <^ Ale ^ ^/X within logarithmic accuracy and the 
function Eq. @, takes the form * = 1 - X'^XlJ2. 

There are two terms in the expression for Y' ^ 
J dQQll i'^'xP + ^Px)^ each of them having a differ- 
ent localization length. The function P determined from 

h>^lddL - >^lcdL + 2Aic5a,. + 2*fcLeo)P = Q (11) 

is given by F = ^Q/(l + ikLco)- This leads to a fast 
decaying contribution to the function Y . The main con- 
tribution at large distances comes from the function 
satisfying the equation 

{n + 2ikL,o)P'x - -XXl.Q . (12) 

Recalling that Q ^ AicAi^, we represent P'^^ in the form 
P'x = XQ\\^f [Xid) and obtain 

(-AL^M. - 2Ai,aAi. + 2ikL,o)f ^ 1 . (13) 

This expression is essentially the same as the one for the 
unitary ensemble, yielding the localization length Leu- 
The pre-exponential of the correlator Y^ is most easily 
estimated for fc = when / = — InAi^. Estimating the 
derivative of Eq. (^) 

Y^^X dAic / dAidlnAid -Xln^X, (14) 

we integrate the result over X and use the fact that the 
correlator Poo(?') at vanishing magnetic field must coin- 
cide with the one of the orthogonal ensemble. This leads 
us finally to Eq. (|). 

It is relevant to note that in our analysis we considered 
the limit when lo Q but the magnetic field H remains 
finite. In principle, this limit differs from the limit H — 0, 
[jj — > considered in Refs. However, for the corre- 

lation function (r), Eq. (|l|), the limit H ^ Q must 
give the same results as those obtained for = 0. 

Although our analysis was performed for very large 
distances r, it allows us to believe that the entire range 
of fields between the orthogonal and unitary ensembles 
is governed by an interplay between the two localiza- 
tion lengths. On the other hand, one should keep in 
mind that Eq. (^) is not exact in the sense that it re- 
mains to be seen how the localization lengths and the 
pre-exponentials change at arbitrary X. 

The localization length Lc enters directly such a physi- 
cal observable as the hopping conductivity . How- 
ever, previous theories as well as the interpretation of 
the experiment of Ref. |^] assumed the presence of only 
one characteristic length. For the hopping conductivity 
in quasi-lD samples, the accounting for the two localiza- 
tion lengths can be performed using Eqs. (^, It is not 
difficult to understand that the activation law, Eq. (|^), 
with Lc — Leo, holds for temperatures (i^irx) ^ < T < 
{viLc) ■ In this regime, the change of the asymptotic 
behavior of wave functions at very large distances is 



not essential. Nevertheless, the far asymptotics is ex- 
tremely important for calculation of the conductivity at 
lower temperatures T < (:/irx) ^. Repeating arguments 
of Refs. [ p^p^ one comes to the conclusion that Eqs. 
(^) can be used also for the description of the latter 
regime provided the temperature Tq is replaced by To/2. 
This means that, experimentally, decreasing temperature 
should lead to a crossover from the activation behavior, 
Eq. (|^) , to another activation behavior with the charac- 
teristic temperature To/2. Measurements at lower tem- 
peratures than those used in Ref. might help to check 
our predictions. 

In conclusion, we have demonstrated that the far tail 
of the density-density correlator in disordered wires at 
arbitrary weak magnetic field decays with the localiza- 
tion length of the unitary ensemble Leu, which is dou- 
ble as large as the length Leo of the main part. This 
means that there is no one-scale crossover between the 
limits H = and H — oo. An arbitrary weak magnetic 
field drastically changes the decay of wave functions: the 
initial decrease with the localization length Leo persists 
only up to the distances rx ''^ —Lc In AT but then it is fol- 
lowed by the decay with the length Leu- An increase of 
the magnetic field leads to a respective decrease of rx so 
that finally, at A ~ 1, wave functions decay everywhere 
with the length Leu- The behavior found can manifest 
itself in the hopping conductivity, which might be the 
simplest way to test our predictions experimentally. 
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